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The field-induced magnetic order in the singlet-ground-state system CsFeCls has been studied by 
measuring magnetization and neutron diffraction. The field dependence of intensity for the neutron 
magnetic reflection has clearly demonstrated that the field-induced ordered phase is described by 
the order parameter {Sx)- A condensate growth of magnons is investigated through the temperature 
dependence of Mz and M^, and this ordering is discussed in the context of a magnon Bose-Einstein 
condensation. Development of the coherent state and the static correlation length has been observed 
in the incommensurate phase in the field region of 5 < i/ < 6 T at 1.8 K. At if > ifc, a satellite 
peak was found in coexistence with the commensurate peak at the phase boundary around 10 T, 
which indicates that the tilt of the c-axis would be less than ~ 0.5° in the whole experiments. 

PACS numbers: 03.75.Hh, 61.12.Ld, 71.70.Ch, 72.10.Di, 75.60.Ej 



I. INTRODUCTION 

Recently, many researchers have revised the field- 
induced magnetic order in the singlet-ground-state sys- 
tem with great interests. Especially, the field-induced 
magnetic order in the dimer spin system have been im- 
pressively investigated in the context of a Bose-Einstein 
condensation (BEG) of magnetic excitations (magnons)^. 
A well-known singlet-ground-state system in the hexag- 
onal compound GsFeCls is constituted by the ground 
state 15'^ = 0) and the first excited states 15'^ = ±1) 
of a fictitious spin 5 = 1 of an Fe^+ ion, each state be- 
ing separated by single- ion anisotropy energy owing 
to a trigonal symmetrical distortion of Gl~ ions along 
the c-axis {c//z). Due to the ferromagnetic exchange in- 
teraction {Ji/kB ~ 4 K) along the c-axis, which is rela- 
tively larger than the antiferromagnetic one in the c-plane 
~ —0.2 K, triangular lattice antiferromagnet), 
the spin system becomes quasi-one dimensional. The ef- 
fect of D is predominant over that of exchange interac- 
tion energies, thus the ground state is nonmagnetic down 
to zero K. However, the magnetic phase transition can 
be induced by applying the external magnetic field along 
the c-axis, then the three dimensional long-range ordered 
(3D-LR0) state appears around the level crossing field 
{He = 7.5 T) of \S^ = 0) and = +1) states at Tn - 2 
K, in the field region of 4 < < 11 T^ The theoretical 
study^ has shown that the phase transition is considered 
to be of the second order and the ordered state is char- 
acterized by the non- vanishing value of {Sx); and also, 
similarity of the present case to the field-induced order 
in a dimer model was pointed out. 

Magnetic excitations and the magnetic phase transi- 



tion in GsFeGls have been extensively investigated ex- 
perimentally by the neutron scattering and NMR exper- 
iments so fari^i^iSiLSiSiiS. The excitation energy decreases 
with lowering temperature, and development of the short 
range order has been found at the wave vector q ^(1/3 
1/3 0) (i^-point)^. Suzuki described these magnetic be- 
haviors by the dynamical correlated effective field ap- 
proximation (DGEFA) theory, which takes into account 
the effect of a fluctuation and correlation via the param- 
eter c»ii. The nuclear magnetic relaxation time Ti of 
^^^Gs has been measured, and the relaxation mechanism 
through the magnetic excitation has been clarified^^^^iiS. 
Comparing the temperature dependence of Ti in the field 
region below and above 4 T, it is suggested that the 
large numbers of magnetic excitations are created and 
a life time of these excitations becomes short in the crit- 
ical region at > 4 T^'^^. Magnetization measurements 
{H // c) have been performed up to 33 T^'^^^ and especially 
the metamagnetism has been paid particular attention. 

As for the 3D-LR0 state, an incommensurate- 
commensurate phase transition has been already stud- 
ied by neutron diffraction up to 5.5 T^{H//c); how- 
ever, the research has been limited in the vicinity of 
the phase boundary, owing to the performance limita- 
tions of cryostats. In this paper, we have focused on the 
magnetic behavior in the 3D-LR0 state, and report the 
neutron diffraction and magnetization studies under the 
magnetic field up to 10 T which completely exceeds He. 



II. EXPERIMENT 

Samples have been prepared in the laboratory at the 
Institute of Advanced Energy (lAE), Kyoto University. 
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The equimolar 99.9% CsCl and FeCl2 powder samples 
were mixed in the ampoule; and after dehydrated us- 
ing HCl gasii, single crystals of a volume of 2~3 cm^ 
large were separated out by the Bridgman method. Clear 
hexagonal cleavage planes were seen, and the color of 
the small crystal (~0.2 cm^) was transparent dark green. 
The mosaic width was about ~ 0.5°, and the hexago- 
nal lattice parameters as used is a=7.237A (the lattice 
parameter of c-axis is not obtained) in the neutron ex- 
periments. 

The magnetization measurements have been per- 
formed using a 9 T SQUID magnetometer at Berlin Tech- 
nical University, which was built on the basis of standard 
construction principles, but supplied with additional fea- 
tures to allow measurements up to 9 T. The measured 
sample was the needle shaped crystal of samplemass 0.61 
mg. The field dependent signal resolution is in the range 
of 5 X 10~^ emu. 

The neutron diffraction experiments were carried out 
on the triple- axis spectrometer El at the Berlin Neu- 
tron Scattering Centre (BENSC), Hahn Meitner Institute 
(HMI), Berlin. The wavelength A=2.4153A was selected 
by Bragg refiection at a pyrolytic graphite monochroma- 
tor for the incoming beam, and the analyser was adjusted 
to measure the elastic part of the scattering. The used 
collimation was 70'-80'-60'-60'. The sample was mounted 
in the new high-field superconducting magnet VMl built 
by Oxford Instruments, which with a maximum field of 
17 T is actually the strongest steady-state magnet avail- 
able for neutron scattering measurement. The prelim- 
inary experiments have been carried out on the triple- 
axis spectrometer PONTA at Institute for Solid State 
Physics (ISSP), Japan Atomic Energy Research Institute 
(JAERI), Tokai. 



III. RESULTS AND DISCUSSION 
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FIG. 1: Magnetization Mz data as a function of external mag- 
netic field (H II c) at 1.8 K. Dissipation of the data at if ~ 9 
T might be the experimental error. 
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A. Magnetization 

First we report the results of the magnetization mea- 
surements. Figure n shows the data at T = 1.8 K, plot- 
ted as a function of field {H He). At such low tempera- 
ture, most of spins should be distributed at the ground 
state; therefore, is nearly equal to zero at low fields. 
Then, Mz shows explicit increase above 4 T in the or- 
dered phase, and is approximately linear up to 9 T. The 
value of magnetization is ~ 1.3 /ie at i^c, which is ap- 
proximately half the magnetic moment expected from 
= 2.54^^. Figure [21 shows the results of the mag- 
netization measurements as a function of temperature at 
5, 7 and 9 T. At 5 T, the broad peak is seen at T 10 
K; this is mainly due to the thermal activation type term 
of exp(— £^i//cb^), where Ei is the gap energy between 
the ground state 15'^ = 0) and the first excited state 
= +1). This broad peak clearly diminishes at 9 T, 
refiecting the change of the ground state from 15'^ = 0) to 
15'^ = +1). Furthermore, the magnetization curves show 



FIG. 2: Magnetization Mz data as a function of temperature 
at 5, 7 and 9 T (i//c). 



bends and the signs of the derivatives change at 5 and 9 
T at T Tn. 

The clear temperature variations of Mz are seen in 
the 3D-LR0 state; however, when each spin is regarded 
as being independent in the molecular field approxi- 
mation (MFA), the magnetization does not depend on 
temperature^. The same feature has been found in the 
dimer spin system TlCuCl^. We expect that the tem- 
perature variations of Mz result from the many body 
effect. Considering that the magnetic order is derived 
from the softening of magentic excitations, the coherent 
state of 16'^ = 0) and 15'^ = +1) will develop following the 
condensation growth of zero energy magnetic excitations. 
Then {Sz) will increase with decreasing temperature in 
the ordered phase by the contribution of 16'^ = +1) com- 
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ponent. 



B. Neutron Diffraction 

Next we show the results of the high field neutron 
diffraction experiments. The previous NMR or neutron 
experiments resulted in failure frequently, since the large 
magnetic field inclines the misalignment of the c-axis fur- 
thermore. In this experiments, the single crystal of a 
volume of ~1 cm^ was attached to the sample holder, 
and combined into the micro-goniomator (Kohzu-Seiki 
CO.); then the obliquity of the c-axis was controlled be- 
low ~ 0.3° at zero field. This configuration was very 
important to observe the incommensurate state^^. The 
former NMR experiments^^ have shown that the ordered 
phase region reduces in the phase diagram and moves to 
the lower fields while increasing the obliquity up to 10°. 

Figure ini shows the dependence of integrated intensity 
on external magnetic field {H // c) for the neutron mag- 
netic refiection at T=1.8 K. We have always confirmed 
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FIG. 4: The dependence of integrated intensity on temper- 
ature for the neutron magnetic reflection at 7.0, 7.5, 8.0 and 
8.5 T {H//c). 
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FIG. 3: The dependence of integrated intensity on external 
magnetic field {H/fc) for the neutron magnetic reflection at 
T=1.8 K. 

the incommensurate state at ~ 5 T, before increas- 
ing the magnetic field to induce the commensurate state. 
The intensity of magnetic (1/3 1/3 0) refiection, which 
is directly related to the amplitude of (Sx)^ drastically 
increased and decreased almost symmetrically with re- 
spect to He. Comparison with the phase diagram.^ clearly 
suggests that the LRO is described by the order param- 
eter (Sx)' It is noted that even a slight tilt of the c- 
axis yielded asymmetry of the scattering intensity with 
respect to i^c^^- Figure^ shows the dependence of in- 
tegrated intensity on temperature for the neutron mag- 
netic refiection at 7.0, 7.5, 8.0 and 8.5 T {H//c). Critical 
point exponents P for the condensation growth have been 
roughly estimated by fitting to the following equation. 



V7 (X (Tn - Tf. 



The values of /3 varied with field, each being 0.30 (7.0 
T), 0.22 (7.5 T), 0.27 (8.0 T) and 0.33 (8.5 T). It should 
be noted, that the value of P at 8.5 T (0.33) is near to 
the value of P (0.35) for a Bose Einstein Condensation of 
the 3D-XY system, in spite of the quasi-one dimensional 
interaction of the system. 

To further investigate the ordered state, we try to esti- 
mate the condensation fraction of magnon. Actually, the 
magnetic structure of the ordered phase has not been de- 
termined yet exactly. From comparison with the isomor- 
phous compound RbFeC^^^, we assume a cone structure 
model for the magnetic structure (120° periodicity along 
the a-axis). Denoting M± and as the antiferromag- 
netic and ferromagnetic ordered moment respectively, the 
magnetic structure factors from each component are de- 
scribed as 



^hkl± I 



^hkl\ 



pU'M-i{l^cos^^)\Ghki\^ (2) 

(3) 



4 



2 ^2 7\^2 ^2 _ _ |2 



^hkl\ 



(1) 



where po is 0.269 x 10~ cm, / is the form factor, 
\Ghki?' = I Yji e^^^^^'^^p is the geometrical structure am- 
plitudes of magnetic ions for the (hkl) refiection, and 
(j) is the angle between the scattering vector and the c- 
axis. The notation of |i^/i/c^±P denotes the fact that the 
magnetic scatterings from the component appear at 
s = Bhki i 'T^^ where B^ki is the reciprocal lattice vector 
and r is the modulation vector for the magnetic struc- 
ture. Because the experiments have been performed at 
T ~ 2 K, we may approximate the Debye- Waller factor 
as 1, then the value of M± can be determined experimen- 
tally comparing the intensities of nuclear (110) refiection 
at zero field (/no; non-magnetic) and magnetic (1/3 1/3 
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0) reflection (/i iq) as following; 
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where = I Ei &i 



iqiio-Vi |2 



is the square of the nu- 



clear structure amplitude factors of all ions. For exam- 
ple, the value of M± is derived as 1.83 and (Sx) as 0.477 
{g^ = 3.84^1) at T = 1.55 K and = 7.5 T. On the 
other hand, the computation of {Sx) in the MFA^ gives 
0.703, when the Hamiltonian parameters are chosed as 
Ji=4.4 K and D=17A K (these parameters are obtained 
from the flt to the experimental values of phase bound- 
ary). The dicrepancy between 0.477 (experiments) and 
0.703 (MFA) is appreciable, but not too large, thus we 
suggest that the MFA theory adequately describes the 
experimental results. 

We notice here again that the obtained values of P and 
M± are rather rough estimations. When the external 
fleld increases io H ^ H^^ there was always a severe 
effect of magnetic torque which increases the tilt of the 
sample. Even though the obliquity was very small, the 
3D-LR0 state as well as condensation of magnon seemed 
to be affected signiflcantly. Moreover, the data points 
are limited in the temperature region of ~ 1 K, where 
the precision of temperature control is ~ 0.1 K. 

So far, we have focused on the symmetrical property 
of {Sx) with respect to i^c, in the commensurate ordered 
state. Meanwhile, spin fluctuation can be asymmetric 
with respect to and such effect might appear on the 
incommensurate state. The origin of the incommensurate 
state has been attributed to the fluctuation effect by P. 
A. Lindgardii and to the dipole-dipole interaction by H. 
Shiba^^. Increasing from zero fleld, the magnetic reflec- 
tion was flrst found around 5 T, at the scattering vector 
q ~ (0.315 0.315 0). Figure [51 shows the fleld variation 
of the scattering vector q and the line width along {hW) 
at T = 1.8 K in the incommensurate phase (5 < < 6 
T). Error bars denote the standard deviation determined 
by fltting to a Gaussian curve. The scattering vector q 
gradually moves to the i^-point with increasing external 
fleld, accompanied by narrowing of the line width and in- 
creasing of the scattering intensity. Development of the 
coherent state and the static correlation length has been 
observed in the incommensurate phase. Under the fleld 
of i7 > i^c, a satellite peak was found in coexistence with 
the commensurate peak at the phase boundary around 10 
T at 1.8 K. Shift of the commensurate-incommensurate 
peak was not obvious. This result can be due to the ef- 
fect of the induced ferromagnetic moment of IS'^ = +1), 
which disturbs the propagation of the magnetic excita- 



tion. On the other hand, if the tilt of the c-axis is more 
than 0.5°, then a satellite peak should have disappeared, 
therefore, it indicates that the tilt of the c-axis would be 
less than ~ 0.5° at such high fleld. 
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FIG. 5: The field variation of the scattering vector q and the 
line width along (hhO) at T = 1.8 K in the incommensurate 
phase (5 < if < 6 T). Error bars denote standard deviations 
determined by fitting to a Gaussian curve. 



IV. CONCLUSIONS 

We have investigated the fleld-induced order in 
CsFeCla in the whole fleld range up to 10 T, by means 
of magnetization and neutron diffraction. The fleld vari- 
ation of the intensity of neutron magnetic reflection has 
clearly demonstrated that the ordered phase is described 
by the order parameter {Sx)- The condensation growth of 
magnon has been investigated through the temperature 
dependence of and M^, and the values of P and the 
condensation fraction of magnons have been experimen- 
tally estimated. Development of the coherent state and 
the static correlation length has been clearly observed in 
the incommensurate phase at 5 < if < 6 T. At if > i^c, a 
satellite peak was found in coexistence with the commen- 
surate peak at the phase boundary around 10 T, which 
indicates that the tilt of the c-axis would be less than 
~ 0.5° in the whole experiments. 
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